Dispersion of tracer particles in a compressible flow 
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The turbulent diflFusion of Lagrangian tracer particles has been studied in a flow on the surface of 
a large tank of water and in computer simulations. The effect of flow compressibility is captured in 
images of particle fields. The velocity field of floating particles has a divergence, whose probability 
density function shows exponential tails. Also studied is the motion of pairs and triplets of particles. 
The mean square separation {A{t)^) is fitted to the scaling form (A(t)^) oc t", and in contrast 
with the Richardson-Kolmogorov prediction, an extended range with a reduced scaling exponent of 
Q = 1.65 ± 0.1 is found. Clustering is also manifest in strongly deformed triangles spanned within 
triplets of tracers. 

PACS numbers: 

Introduction. — Many processes in nature depend on the efRcient diffusion in a turbulent medium such as spreading 
of drifters on the ocean surface P, I3 or of aerosol particles that can act as nuclei for the formation of rain droplets in 
the atmosphere If the flow is compressible (or the particles have inertia), the following can occur: stirring an 

initially uniform distribution of particles can induce them to coagulate, as can be seen in Fig. 1 (left). This photo was 
taken a fraction of a second after floating particles were uniformly dispersed on the surface of a large tank of water 
(lateral dimensions 1 m x 1 m) that was continuously sustained in a turbulent state by bulk driving. The particles 
appear in white in this figure. After a few seconds they accumulate in narrow ridges with voids in between, forming 
a lace-like network. 

The origin of the coagulation is apparent in the right panel of Fig. 1. The incompressible fluid below the surface S 
is shown to be moving upward near point U and downward near point D. The floaters cannot follow the bulk motion; 
rather they accumulate at points like D and flee points like U. The floating particles move in the x — y plane, which will 
be taken at z = 0. They simply sample the horizontal components of the flow velocity, Vx{x,y,0,t) and Vy{x,y,0,t), 
of the incompressible fluid below them. Since water is incompressible, the divergence of the velocity field is zero at 
all points on the surface and below it, so that dvx(x,y,z,t)/dx-\-dVy{x,y,z,t)/dy = —dvz{x,y,z,t)/dz, which is 
non-zero and, in fact, large. Thus the floaters form a two-dimensional compressible subsystem, even though the flow 
velocity is negligible compared to the speed of sound and their behaviour is unrelated to two-dimensional turbulence 
i 




FIG. 1: Left: Image of a cloud of particles on the surface of a turbulent body of water. The particles which are 10 /^m in 
diameter appear as white and were initially dispersed uniformly over the surface. This image was captured 100 ms later. The 
scale below the figure is in cm. Right: Sketch illustrating the origin of the clustering effect at the surface S. 



2 



Recently, compressible turbulence has come under theoretical scrutiny (see [a for a review). Modelling the flow as 
Gaussian and delta-correlated in time, one can show for both spatially smooth Q and spatially rough flows Q, that a 
transition from a regime of weak to strong compressibility occurs at some threshold value of the dimensionless ratio, 
defined as 



C EE {{dvjdx,)^)/{idvjdx,)^) , (1) 

with i,j — 1,2. Clearly, for any C ^ tracers will cluster, but the eventual particle distribution will depend on 
whether C is above or below the threshold. Assumptions contained in those models make it inapplicable to the present 
experimental situation, where the flow is non-Gaussian and has a finite time correlation. Because the idealizations in 
those models are not realized in the experiment, we do not know if such a transition from weak to strong compressibility 
exists in our experiment. On the other hand, a computer simulation (DNS) that roughly matches the flow conditions 
of this experiment (incompressibility of the three-dimensional velocity fleld at all z < 0, flnite correlation time,^nd a 
free-slip boundary condition at the flat surface) are in remarkably good agreement with the first experiments [3, [l3 ■ 
Herein, we summarize and extend recent experiments using freshly cleaned water surfaces [T]| . 

Our study focuses on two points. The first one is the relative motion of the floaters, with special emphasis on so-called 
Richardson diffusion which refers to the time dependence of the mean square particle separation (A(t)^) ^ t" 
with a = 3. Scaling close to the Richardson prediction was found in experiments [iSj as well as in high-resolution 
numerical simulations in the inverse cascade of two-dimensional incompressible turbulence. Intermittency effects 
if present seem to modify this exponent only slightly fl^. Dimensional arguments would seem to be inapplicable 
for the floaters, since they can exchange kinetic energy and mean square vorticity with the fluid below them. Since 
a theory for their motion has yet to be developed, one must again turn to computer simulations, which have been 
recently carried out p^ . 

Secondly, first studies the competing effects of coherent shear and random motion in the surface flow are presented 
by analyzing the evolution of shapes formed by three particles tracked simultaneously 0, . By neglecting the 
center of mass, the relative position of three points can be conveniently expressed by the following two vectors, 
ai — (x2 — xi)/-\/2 and a2 = (2x3 — xi — X2)/-\/6 0|. The distortion of a triangle is measured by the ratio /2 = g^/i?^ 
which relates its radius of gyration, R = \/ a\ + a|, to the smaller half-axis, 527 of the smallest ellipse that covers the 
tracer triangle. From il^ we get for two dimensions 
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We will present measurements of the evolution of (J-i) as a function of time and propose a qualitative explanation 
based on the compressibility of the flow. 

Experimental setup. — The square tank in which the measurements were made was fllled to a depth Z, which was 
typically 30 cm, and the turbulence was generated by a variable speed 8 hp pump that circulated the water through 
a series of pipes capped with rotating jets located well below the surface. The turbulent intensity was limited to keep 
the amplitude of the surface waves below 1 mm. Images from a fast overhead camera, which typically operated at 
300 frames/s, were captured in a computer and analyzed with in-house software. A 5 W diode pumped solid state 
laser was used to produce a sheet of light illuminating only the surface. 

Measurements of the velocity well below the water surface yielded a second order longitudinal structure function 
D2{r) (X r'^-^^^o.od ^j^^ interval 1.2 mm < r < 3.5cm. This is close to the Kolmogorov value of 2/3. On the 
surface, the outer scale of the turbulence rg was 3.2 cm. Essential to the success of the experiment was the continuous 
vacuuming of the surface by a small auxiliary pump connected to a vertical pipe that comes up to the surface. Without 
this cleaning, which was already reported in jili] . the particles will begin to interact as their density becomes large, 
destroying the clustering seen in Fig. 1. The floating particles used in velocity measurements on the surface were 
hollow glass spheres with radius a = 200 fim and highly buoyant particles of specific gravity 0.25 with a =50 fiui. The 
Stokes time, ts = 2a^/(9z^), was roughly 1 /xs, so that tracers easily follow the fiow beneath them {i/ is the kinematic 
viscosity). Other particles, such as mushrooms spores and talc were also used successfully. 

Surface flow divergence. — By tracking the motion of the small surface particles over a short interval of time, one 
can map out the velocity Vx{x,y,0,t),Vy{x,y,0,t) of the floaters and extract the local surface divergence, V2 • v — 
dvx{x,y,t) / dx + dvy(x,y,t)/dy, a quantity that will fluctuate in time and space. The two-dimensional divergence 
V2 • V is therefore a random variable, whose probability density function (PDF) Pva v is plotted in Fig. 2. Those 
measurements were made at Taylor microscale Reynolds numbers Re of 100 and 140 (left panel). The relatively 
large size of the seed particles made it impossible to directly measure velocity differences on a scale r less then the 
Kolmogorov dissipation scale, 77, which is 0.1 mm. Therefore, the quoted Taylor microscale was determined from 
gradients measured over a few millimeters. 
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FIG. 2: Probability density function (PDF) of the surface flow divergence field with a — ((V ■ v)^)^''^. Left panel shows the 
experimental findings for two Taylor microscale Reynolds numbers of Re — 100 and 140, respectively. Right panel shows the 
simulation results for three different values. Runs at Re — 106 and 157 are for driving by a large-scale shear. Data from [l^ 
were obtained with isotropic volume forcing at a Taylor Reynolds number of Re — 145. 



New experimental data were extracted from an ensemble of measurements made at many instants of time and many 
points {x, y) in the surface. Note that the mean value of V2-v is close to zero and that this function has exponential 
tails out to almost four decades in its standard deviation, CTVa-v- Computer simulations of P^^ v for isotropic forcing 
[l^ as well as large-scale shear likewise show this exponential behaviour (see the right panel of Fig. 2). Exponential 
behaviour for the PDF can be expected, because the divergence is composed of velocity derivatives which are known 
to yield exponential tails in two- and three-dimensions. Note that the slope of the decay of the tails is effected by the 
details of volume forcing and seems thus to be non-universal. In both the computer simulations and the laboratory 
experiments, it was found that C « 0.5 regardless of the intensity or geometry of the forcing. 

Pair dispersion. — An important characteristic of turbulence is the ability of velocity fluctuations to disperse par- 
ticles. Figure 3 displays relative diffusion (A^(i)) for both laboratory experiments (circles) and computer simulations 
(triangles). One expects that a =3 for an incompressible turbulent flow and times larger than the Lagrangian decor- 
relation time which is ~ 10 Kolmogorov times, t^. Instead, both the simulations and the experiments on the present 
flow show a crossover from a ~ 2 to an exponent in the range 1.65 to 1.8. Although the simulations display a small 
range where the slope is approximately 3.2, the most striking feature of the pair dispersion measurements is a reduced 
value of the scaling exponent a in the inertial-range for both measurements and the simulations. 

The reduced a reflects a correlated motion of particles within a pair: (i) pair particles remain close together while 
moving in a big cluster or pair trajectories merge again after longer transient phases of some separation, i.e. their 
contribution to the growth of (A^(i)) with time is small, (ii) on the other hand, decorrelation due to breaking of 
clusters by newly upwelling fluid or simply large scale separation does appear. This break-up introduces a randomizing 
element and causes a diffusive Brownian scaling. While the latter mechanism can be found in an incompressible system 
as well, the first scenario, (i), might be related to the compressibility of the advecting flow. Note that one has to be 
careful with the interpretation of our results. For large compressibilities, trapping of trajectories can be predicted 
for the Kraichnan flow analytically [sj. Recently, Sokolov performed numerical particle dispersion experiments in a 
synthetic incompressible flow where the exponent P of the spatial scaling of the distance dependent correlation time 
of the quasi-Lagrangian pair velocity difference could be varied, t(A) ~ A'^ ^M- F'or /3 increasing from to 2/3 
the exponent a of the pair dispersion increased from 3/2 to 3. It is still unclear how the compressibility affects this 
correlation time scaling although flrst attempts have been made |20l |. The present resolution of the data prohibits a 
detailed analysis of this issue for a Navier-Stokes flow. 

Experimentally, is determined from t^~^ = (dvx/dx) measured at rj. The value of Tjj was extrapolated from 
measurements at larger scales and was found to have a maximum value of 1/80 sec. The measurements of the 
surface velocity field, Vx{x,y,0,t),Vy{x,y,0,t), enable determination of the mean square pair separation A^(Ao,i), 
where Aq is the pair separation at i = 0. It was not possible to accurately track the motion of true floaters over 
large distances because of the slight vertical motion of the surface and the time-varying illumination it produces. 
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FIG. 3: Pair dispersion {/S^{t))/'if vs. t/r,, for the surface flow experiment (open circles). For a better comparison we added 
the result of the numerical experiment (open triangles) p^ . 



Instead, the measured velocity fields were used to propagate the motion of simulated particles where we checked the 
robustness of the pair diffusion measurements with respect to variations of the interpolation parameters. Snapshots 
of the simulated particle fields verify that the structures produced are qualitatively the same as for the real particles 
seen in Fig. 1. The limited resolution of the data prevented an analysis of higher order moments of the particle density 
which is another appropriate measure to characterize the clustering (for more details see |,16:j). 

Shape distortion of tracer triplets. — In order to further probe the effects of compressibility on flow structure, we 
investigated the distortion of triangular configurations of floaters . The experimental data for I2 are in the left upper 
box of Fig. 4 and the simulations are shown in the box to the right. In both cases the particles were initially placed 
to form equilateral triangles, for which I2 = 1/2 (i.e. ai _L a.2 and |ai| = |a2|). The three curves are for triangles of 
different initial side length. 

Pumir et al. |17| argue that coherent shear is responsible for flattening initially equilateral triangles, causing 
a decrease in (/2)- After a sufficiently long time, random small-scale fluctuations tend to drive {I2) towards the 
Gaussian equilibrium value of I2G = (1 ~ the two-dimensional incompressible case, the experiments show 

that the Gaussian limit is achieved in roughly 20 t/r^ [l^. It was asserted in Ref. Tt'I that coherent shear should 
dominate in the viscous subrange and give way to the randomizing effects of turbulent mixing at larger scales. 

The measurements and the simulations in Fig. 4 show a strong dip below the asymptotic value of I2G and an extended 
minimum that exists for triangles initially having a radius of gyration significantly larger than the dissipative scale. 
This behaviour can be accounted for by strong clustering that is apparent in Fig. 1. The triplets are quickly compressed 
into quasi-colinear configurations, making ai almost parallel to a2 which forces /2, shown in Fig. 4, towards zero (see 
Eq. 0). In incompressible turbulence, this effect is absent, and I2 very soon attains its random Gaussian value l2G- 
The laboratory data attain a value of I2 that is slightly larger than the Gaussian value. The triangles in the computer 
simulations show an even stronger tendency to remain flat and appear to approach an equilibrium value smaller than 
l2G- Self avoidance of the finite-size particles and possibly the particle tracking resolution used in the laboratory 
experiments force a limit on the minimum attainable value of {I2) and presumably cause its later growth. In the 
numerical study (point) particles can become arbitrarily close and thus {I2) can remain at a small value well below 
Gaussian. 

The above results from the computer simulations suggest that the compressibility is acting to flatten triangles, 
even when R is much greater than 77, contrary to incompressible flows [T^ . A dominant number of triangles remain 
distorted for very long times. This is illustrated in the two lower panels of Fig. 4, which are plots of the joint PDF 
of the normalized vector magnitudes |ai|/i? and |a2|/i?. After a short time (lower left panel), the triangles are only 
slightly distorted, which is reflected by weak scattering about the equilateral value, |ai|/i? — |a2|/i?. The right 
panel shows that the triangles have become strongly distorted at a sufficiently later time. Here we observe two peaks 
(fli ^ 02 and vice versa). It is notable that the distribution function is confined within two sharply defined peaks 
rather than spread widely through the |ai|/i?- |a2|/i? plane. Laboratory data pertaining to these distributions is 
not reported here. However the analogous PDF's have also been measured. The finite size of the particles blocks 
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FIG. 4: Upper row: Temporal evolution of (/g). Left panel is for experiment at a Taylor Reynolds number of Re = 100. Right 
panel is for numerical simulation at Re = 145. The initial side lengths of triangles are given in the legend. The Gaussian 
asymptotic value of the two-dimensional incompressible case of I2G = (1 — '^1^)1'^ is indicated as a dotted line. Lower row: 
Probability density function p{ai/R,a2/R) from the simulations for two points in time and the initial triangle side length of 
1377. Both times arc indicated by filled squares in the upper right panel. The support of the PDF is the unit circle. The initial 
equilateral case would correspond to a delta-peak at 45 degrees. This distribution spreads toward the cases where ai S> 02 and 
vice versa. 



observation of the peak at |ai|/_R = and |a2|/i?=l, apparent in the simulations, in the lower right panel of Fig. 4. 

In summary, the relative motion of passive tracer particles in our compressible flow is very different from tracer 
motion in incompressible fluid turbulence. The turbulence creates locally convergent and divergent regions of par- 
ticle density and inhibits Kolmogorov-Richardson pair dispersion. Likewise, the evolution of geometrical structures, 
revealed in the motion of particle triplets, strongly differs from incompressible fluids. 
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